AC field induced quantum rectification effect in tunnel junctions 
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We study the appearance of directed current in tunnel junctions [quantum ratchet effect), in the 
presence of an external ac field fit). The current is established in a one-dimensional discrete in- 
homogeneous "tight- binding model". By making use of a symmetry analysis we predict the right 
choice of f{t) and obtain the directed current as a difference between electron transmission coeffi- 



cients in opposite directions, AT = T 



T . Numerical simulations confirm the predictions 



of the symmetry analysis and moreover, show that the directed current can be drastically increased 
by a proper choice of frequency and amplitudes of the ac field f{t). 



PACS numbers: 05.60.Gg, 73. 23. Ad, 73.40.Gk 

A great attention has been devoted to theoretical and 
experimental studies of a transport rectification in vari- 
ous physical, chemical and biological systems In 
particular this peculiar effect appears in the form of a di- 
rected current as a system is exposed to a time-dependent 
ac force with zero mean. Such a directed motion (cur- 
rent) has been observed in molecular motors 0, 13, in 
Josephson junction coupled systems 0, 0|, and in sys- 
tems of cold Rb atoms in laser fields p|, just to name a 
few. 

In most studies the appearance of directed current has 
been analyzed by making use of the model of a particle 
moving in thepresence of both space and time periodic 
potentials S B 13 • The effects of dissipation and in- 
teraction with a heat bath have been taken into account, 
and in many cases the decrease of dissipation as a system 
approaches the Hamiltonian limit, leads to a substantial 
increase of the directed current 9] . By making use of a 
symmetry analysis of the dynamic equations of motion 
it has been shown that the directed current occurs as 
a space-periodic potential is asymmetric (ratchet poten- 
tial) and/or all important symmetries of an ac force, i.e. 
time-reversal and shift symmetry, are broken 0, |^ ■ No- 
tice here, that the latter case can be easier realized in 
experiments. It has been established that more complex 
systems described by generic Bolzman transport equation 
[lo|, Fokker-Planck equation 0,IEE|i etc. also display 
a directed current. Moreover, the directed energy trans- 
port has been obtained in the case of interacting many 
particle systems described by nonlinear partial differen- 
tial equations Q ■ 

Most previous studies analyzing both the dissipative 
and Hamiltonian limits, were heavily based on a classi- 
cal nonlinear regime of particle motion. A next step 
is to obtain the directed transport in systems displaying 
quantum-mechanical behavior. The quantum regime of 
directed transport for a few particular systems has been 



studied in Refs. 113 . Iisl H^ . It was shown that a 
dissipationless quantum system in a single energy band 
tight-binding approximation does not support directed 
current 14], and although a quantum Brownian particle 
can display a directed motion and even numerous cur- 
rent reversal in the presence of an asymmetric ratchet 
potential, one still needs to take into account the dissi- 
pative effects Notice that a quantum ratchet effect 
has been observed in Josephson junction arrays as 
the directed motion of a Josephson vortex in a ratchet 
potential with two interacting energy bands. 

On the other hand there is a class of intrinsically quan- 
tum phenomena involving a particle tunneling through a 
potential barrier where effects connected with the dissi- 
pation in tunnel region can be diminished. The particu- 
lar examples of systems displaying these phenomena are 
the ballistic regime in mesoscopic devices, various tunnel 
junctions, tunneling through quantum (molecular) wires, 
etc. Recently the directed current has been obtained for 
electrons that tunnel through the barrier containing the 
asymmetrically distributed energy levels in the presence 
of harmonic ac drive [l^ llVj. However, in this case the 
dissipation caused by the coupling to the leads is crucial 
in order to obtain a nonzero directed current. Thus, it is 
interesting to obtain general conditions for the appear- 
ance of directed transport in tunnel (dissipationless in the 
tunnel region) systems that are subject to an externally 
applied ac drive. 

In this Letter we consider electron tunneling through 
a one-dimensional potential barrier U{x) in the presence 
of an externally applied ac force f{t). Moreover, we will 
restrict ourselves to the elastic scattering case, i. e. the 
energy of incoming electron coincides with the energy of 
transmitted electrons. By making use of symmetry argu- 
ments we obtain the conditions for the directed current 
in such a system and verify these arguments by direct 
numerical simulations. We show how the directed cur- 
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rent can be resonantly enhanced by a proper choice of 
the frequency and ampUtudes of ac drive. Notice here, 
that the obtained current can be considered as a particu- 
lar^exajnpl e of the widely discussed quantum pump effect 

In the absence of dissipation the propagation of elec- 
trons through the potential barrier U{x) (U{x)x^±oo — > 
0) is described by the one-dimensional Schrodinger equa- 
tion. The externally applied ac force f{t) induces a time- 
dependent potential that is odd in space, i.e. —xfit). 
The crucial characteristics of electron propagation is the 
transmission coefficient T{E) which in the elastic case 
depends on the energy of electron E. Moreover, we will 
distinguish the electron transmission coefficients in the 
opposite directions: T^^{E) as the electron moves from 
left to right and T^^{E) as the electron moves from 
right to left. Thus, the directed current in such a sys- 
tem appears if a symmetry is broken and the difference 
AT{E) = T^^{E)-T'^^{E) is not zero. The symmetry 
properties of the transmission coefScient T{E) can be an- 
alyzed by making use of the Green function method jl^l 
and es pec ially its representation in the form of path in- 
tegral (23. Indeed, the transmission coefficient T{E) is 
written in the form: 



TiE) 



oc 



— t £/r 

dte G[xi,X2;t) 



(1) 



where the coordinates xi and X2 are located on the dif- 
ferent sides of the barrier. Thus the symmetry of trans- 
mission coefficients T'"^'^^^\E) and correspondingly the 
absence of directed current is determined by a symme- 
try of the Green function G{xi, X2', t) with respect to the 
permutation of xi and X2. In order to obtain the condi- 
tion allowing a nonzero directed current we represent the 
Green function in the form of the path integral [2^: 



G{xi,X2]t) 



D[u(T)]e" 



D[u(t)]e'' 



iS{u{t): XI, X2) 



-{'^r-U{u{T))+u{r)f{T) 



(2) 

Here, m is the electron effective mass, and all paths u(t) 
start from the coordinate xi at t = — 00 and pass the 
point with the coordinate X2 si t = t. A most important 
symmetry transformation is time reversal. We notice 
that time reversal t — > — r changes (permutes) the co- 
ordinates xi and X2, and the Green function G(xi, X2', t) 
is invariant with respect to this permutation and corre- 
spondingly the directed current is absent, provided f{t) 
is symmetric [23): 



fit) = /(-i) 



(3) 



A simple consequence of such a time-reversal symmetry 
is that the directed current vanishes for an arbitrary po- 
tential U {x) in the absence of an ac force. 



Next, we consider another symmetry operation: re- 
flection in space. In this case the symmetry operation 
X —>■ —X, and correspondingly u — > —u in the integral of 
Eq. (O, changes the coordinates xi and X2 but leaves the 
Green function G{xi, X2',t) invariant, provided the poten- 
tial U{x) is symmetric and the ac force is shifl symmetric: 



U{x) = Ui-x); fit) = -f{t+-) 

UJ 



(4) 



where oj is the period of ac force. 

It is interesting to mention that the conditions Q and 
Q of the absence of directed transport obtained for the 
quantum-mechanical system, coincide with ones obtained 
for classical dynamical systems in Hamiltonian limit 0, 
Thus, a simplest choice to obtain the directed 
current is to use the symmetric potential U{x) and the 
ac force f{t) allowing to break the shift symmetry and 
time-reversal symmetry 



f{t) = /i sinwt + /2 sin(2a;i + (^). 



(5) 



Note that only for the phase shift (p = ± 7r/2 the ac 
force f{t) possesses the time-reversal symmetry, and the 
directed current has to vanish. 

To verify the above presented analysis and also to 
study the dependence of the directed current on the fre- 
quency and amplitudes of the ac drive, we consider a par- 
ticular example of a tunnel junction, namely a large and 
wide barrier between two narrow-band electrodes (Fig. 
1). The ac force of the type ||SJ) is applied to the barrier 
region. In a simplest case this system can be described 
by a standard tight-binding model in the presence of a 
time-dependent scattering potential: 

ifltpn = Co(2'0„-^„+i -'0„_i)-|-C/„'0„-/„(t)l/'n . (6) 

The potential C/„ and /«(<) have the form: 



Un 



Ut) 



Uo , -2 < n < 2 
, n < -2; n > 2. 



-nf{t) , -2 < n < 2 

, 71 < -2; n > 2. 



(7) 



(8) 



For simplicity the hopping between the sites determined 
by Co, is assumed to be the same in the electrodes and 
in the barrier region. Note here, that this model can be 
directly used also for the analysis of electron propagation 
through the periodic chain of impurities with symmetri- 
cal (see, Fig. 1) or asymmetrical distribution of energy 
levels [Illia. 

We use Co = 0.1 to ensure that the case of elastic scat- 
tering (a one channel scattering) was realized. The height 
of the barrier was fixed Uq ~ 1 and thus, in the ab- 
sence of ac drive the transmission coefficient was negligi- 
bly small ( ~ 10~^^). In order to resonantly enhance the 
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FIG. 1: The schematic view of a system (tunnel junction) 
allowing to obtain a quantum-mechanical directed current. 
Dashed lines show the boundary of a tunnel barrier region. 
The tunnel junction is modelled by a chain of different atoms 
(black circles). The electron hopping between atoms is the 
same in electrodes and in the barrier region, but the energy 
levels of atoms in the barrier region are shifted and form a 
tunnel barrier. The ac force is applied to the barrier region 
only. 



transmission coefficient and the directed current the fre- 
quency = 0.5 was used. For this particular value of a; 
both harmonics in |SJ) give a contribution to the resonant 
enhancement. Such an enhancement of the transmission 
coefficient is due to a well known photon assisted tunnel- 
ing. As the ac drive is applied, the electron can absorb 
one or two photons, then propagate through the barrier 
region with the enhanced energy E = E + 2huj > Ua 
and emit the phonons at the end of this region. The 
transmission coefficient T{E) for such a process can be 
large and moreover, can be different for electrons propa- 
gating in different directions. 

To calculate the transmission coefficient wc turn to a 
numerical investigation. We developed a scheme that 
generalizes the one given in E^l, and allows to compute 
the transmission coefficient for any time-periodic scatter- 
ing potential. The details of this method can be found in 
|25| . and this scheme has been used in order to calculate 
the phonon scattering by time-periodic and localized in 
space breathers |25j . 

The dependence of the computed transmission coeffi- 
cient T^^{q) on the wave vector q of electron propagating 
from left to right is shown in Fig. 2. Notice here that 
in our model the wave vector q is determined by a stan- 
dard relationship: E = 4co sin^(g/2). As expected the 
transmission coefficient drastically increases and more- 
over, displays the effect of resonant tunneling for partic- 
ular values of the energy E' of a propagating electron. 
For these values of energy the transmission coefficient is 
close to one. Notice that the resonant transmission ob- 
tained for the electron wavevector g — tt/ 2 has a simple 
physical origin. Indeed, for this particular energy the 
resonant condition qL = 27r, where L is the length of 
the barrier region, is satisfied and the transmission has 
a maximum. As the amplitudes of ac force /i^2 or the 
phase shift Lp vary additional resonances appear in the 
dependence T'"^{q) (see Fig. 2). These resonances in 
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FIG. 2: The wave vector q dependence of the transmission 
coefficient T^"-{q): one harmonic case /i = 0.9 and /2 = 
(soUd line), and two harmonic case where the amplitudes 
/i = 0.9 and /2 = 1 were used (dashed line). The value 
If = TT /2 was used. 
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FIG. 3: The dependence of AT on the phase shift ip. 
The curves are shown for a particular value of wave vector 
q — n/2 as the resonant transmission occurs. Solid line: 
/i = 2.36 and ji = — 0.03; dashed line: /i = 2.36 and 
/2 = — 0.003. Note here that the latter case is scaled by 
a factor of 10, so the real values of AT is 10 times less than 
they appear in the plot. 

T^^{q) are due to the resonant propagation of an elec- 
tron with enhanced energy E — E + 2huj through the ac 
force dependent energy levels of a barrier region. 

Moreover, we observed a strong dependence with a sign 
change of the difference between transmission coefficients 
AT{E) = T^^(E) - T^^(E) on the phase shift ip 
(see Fig. 3). As it was expected from generic symme- 
try analysis, AT{E) is practically zero as the phase shift 
(fi — ± tt/2. The value of AT{E) strongly changes with 
the amplitudes of ac drive, and the maximum value of 
AT{E) ~ 0.043 (around 7.5% from the value of T{E)) 
has been found for particular values of /i^2- 
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Notice, here that an increase of the spatial extension 
of the ac drive " spot" can lead to a substantial decrease 
of both transmission coefficient and directed current. It 
is due to a specific of our model where a small value of 
coupling Co has been chosen. 

The physical origin of the considered rectification effect 
is the quantum-mechanical interference between various 
paths of an electron with different absorption (emission) 
of photons. For this particular case of ac induced reso- 
nant propagation in tunnel junctions, this interference is 
especially transparent. As an electron propagates from 
left to right it absorbs one or two photons on the left 
side of a junction, and after the propagation, emits the 
photons on the right side of the junction. The ampli- 
tude of this process is A{-f^ + if2e'^){-fl + ihe^'^f): 
where A is the amplitude of electron transmission with 
an enhanced energy. Similarly, for the electron propagat- 
ing from right to left we obtain the amplitude A(—f1 — 
if2e-'^'^){—fi — if2e~'^'^), and correspondingly, the differ- 
ence in the transmission coefficients occurs, which is pro- 
portional to cos{(p). 

In order to observe this resonant rectification effect 
the ac drive of the frequency that is a half of the barrier 
height has to be used. Moreover, the second harmonic 
also has to be present. Thus, it can be obtained with the 
standard infrared laser sources if the barrier height Uq is 
around Q.leV. Our results can be used not only for tun- 
nel junctions but also for the wave propagation through 
various transmission lines, e. g. Josephson junctions ar- 
rays, in the presence of both a localized static scattering 
potential Uq and an externally applied mixed ac drive of 
frequencies lo = Uq/{2K) and 2w. 
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